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Abstract. - It is shown that magnetic bubble films behaviour can be described by using a 
2 D super-Ising hamiltonian. Calculated hysteresis curves and magnetic domain patterns are 
successfully compared with experimental results taken in literature. The reciprocal problem of 
finding paramaters of the super-Ising model to reproduce computed or experimental magnetic 
domain pictures is solved by using a perceptron neural network. 



Thin films with perpendicular magnetization show widespread technological applications, 
particularly as high density storage devices. Thus, our understanding, at a fundamental level, 
of the way the magnetic domains are organized in such films is extremely important. This 
letter is focussed on films having cylindrical shape bubble and labyrinthine magnetic domains. 
These films have been extensively studied in the seventies as they were considered as a possible 
support for high density magnetic recording, essentially through the possibility of manipulat- 
ing bubbles [I]. There is nowadays a renewed interest for bubble and stripe domains [2,3], 
because of the appearance of new techniques, like Magnetic Force Microscopy (MFM) or 
Secondary Electron Microscopy with Polarization Analysis (SEMPA), showing these domain 
geometries at the nanometer scale. The theoretical models used up to now to describe these 
domains were developped 30 years ago [4-6] , and consist of a comparison of the energies of do- 
mains having different ideal geometries. Such an approach is the simplest one can imagine but 
provides a limited understanding of the mechanism of domain formation. More elaborated de- 
scriptions based on linearized micromagnetic equations [7] or Ginzburg-Landau formalism [8] 
have been proposed. These approaches describe qualitatively the transition from bubble to 
stripe structure, but they do not allow quantitative comparison with experimental results. 

It will be explained in a first part that bubble films behave as an Ising system with long- 
range interaction (super-Ising) between bubbles considered as giant spins. The influence of 
the interaction potential parameters (range and shape) on the domain geometry and on the 
magnetization curve will be discussed. The comparison between simulations and experimental 
results will help to determine the relevant parameters. In a second part, we solve the inverse 
problem i.e the determination of the interaction potential from a simulated or experimental 
magnetic map. The treatment proposed is a perceptron type neural network (NNT) searching 
for correct values of the potential using an error-correction learning rule [13]. To our knowl- 
edge, it is the first time that neural network techniques are proposed to analyse magnetic 
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domain pattern geometries. 

The magnetic behaviour of bubble films can be understood starting from a magnetic mono- 
domain situation. This state becomes unstable if an initially large perpendicular magnetic 
field is decreased to reach a critical value Bn (nucleation field). In fact, for such a field, 
a situation with lower energy is obtained if a bubble with a diameter is formed (critical 
bubble). The competition between the wall and magnetostatic energies determines do. The 
magnetostatic energy favours the formation of domains and varies, in first approximation, 
with the volume (Vb oc dg) of the bubble. The wall energy, that does not favor the presence 
of domains, is proportional to the bubble's surface (5*0 oc do). The wall thickness is assumed 
to be negligible [1-6]. The characteristic length I = cr^/47rM^ refiects the balance between 
the wall energy and the magnetostatic energy, cr^, is the wall energy density and Mg is the 
saturation magnetization of the material. Typically, for cobalt or iron-palladium alloy, I is 
around 10 nm. This length allows us to treat a non-dimensionnal problem by choosing 
AttMs and AirMghdQ as respectively length, field and energy units. The function do (h) has 
been tabulated [1,5] and can be approached by do = 0.15ft. + 3.86 for h < 50 ; the tilded 
quantities are non-dimensionnal variables. 

We have considered up to now an isolated bubble. The particular case of two down mag- 
netized critical bubbles in an up magnetized film (see fig. 1) is the following. Either the two 
bubbles are not in contact and they only interact via the dipolar field, the wall energy is then 
2 X aw So, or the bubbles are in contact and the wall energy becomes aw {"2 x 5*0 — ^i). Si 
being the contact surface between the two bubbles. These two cases can be summarized in 
a single model on a square lattice in which each cell i, of size do x do is given the value of 
r]i = ±1, corresponding to up and down magnetization. Then the interaction energy between 
two cells reads: 



rpi'. 



E''''^{\^^3\) mVj (1) 



where the first term is the wall energy and the second term the dipole energy. The hamiltonian 
fully describing the interaction between cells is given by: 

with Vij = 5j^i±i/2do — E''-^p {\i — j\) if i ^ j and Vi^i — 0. The second term accounts for 

the Zeeman coupling between each bubble's momentum and the applied magnetic field B. 
The hamiltonian ^ is formally identical to a super-Ising one. Similar analysis, using the 
competition between wall and dipolar energies, were previously presented in [9, 10] but it is 
the first time that bubble magnetic films are shown as 2 D super-Ising systems. It can be 
noticed that a similar hamiltonian could be written to describe perpendicularly magnetized 
nanoscale dot arrays [11]. 

If the dipolar field of a square cell is approximated by the one of a cylindar of radius 
Ro = do/2, E'^'P reads: 

^..P i\^-J\)^^ j^^ de (:^) ' Jo (. V - J\) (l - e--/«") (3) 



where Jo and Ji are Bessel functions of order respectively and 1. The dipole term is short 
ranged even if its range increases with the film thickness (see fig. 1). For simplification of the 
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Fig. 1 - Interaction dipolar energy between two bubbles of same sign as a function of their separation 
distance R for different film thicknesses. The inset shows two down cylindar shape bubbles in an up 
magnetized film. 

forthcoming discussion, the potential V is approximated by: 

^^'~\ y_ if i?+<|z-j|<i?_ ^ ^ 

where V+ is the sum of the dipolar and the wall energy between first neighbours, is the 
dipolar energy, R^^ the distance to the first neighbours and i?_ is the range of the potential. 
For the cobalt films studied in [2], with thickness h — 50 and do — 10, V+ and V_ can be 
respectively fixed to 0.02 and -0.02, = 4 and i?+ — ^/2. This choice of R+ is equivalent 
to consider also the second neighbours of a cell as nearest neighbours. This allows the model 
to be more isotropic and attenuates the lattice effects. 

With Ai = V (k) r]i+k + B, the average momentum of a cell i, for a given environment 
enwi, is in the canonical ensemble: 

{Vr)env, = t^nh ^J^" \j (5) 

Since the ratio V+/kB ~ 10^ K, in a reasonnable T range, one can let T go to 0+ in 101 and 
let the lattice follow the rule: 

^^^0 (6) 

[ random choice if A; — 

The case A^ = corresponds to rji equals ±1 being equiprobable. This is crucial since 
decreasing B from the monodomain state (all rji = +1), the magnetization starts to flip for 
Bpf = — ^i;V{k), i.e. Ai — 0. Without dipolar field, the magnetization sharply flips for 
Bn < with no domain nucleation (in this case i?jv is similar to a coercitive field). With 
dipolar field, Bjv is possibly positive and, in this case, the magnetization begins to decrease 
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Fig. 2 - Simulated magnetization curve. The computed patterns (100 x 100 cells) show the domain 
evolution at different points of the hysteresis curve when the field is decreased, starting from the 
saturated up state. Black and white parts correspond to up and down domains. 

Fig. 3 - Computed patterns in the demagnetized state for different potentiel values, R+ — \/2 and 
R- = 4. Black and white parts correspond respectively to up and down domains. 



by creating domains: Sjv is the nucleation field. 

The hysteresis curve shovifn in fig. 2 is calculated for the parameters defined above. The 
similarity with the experimental results from [2] is striking. Our computed value of Bjy ~ 10 
kOe (Ms = 1480 erg.cm"^ for cobalt) is in good agreement with regard to the experimental 
value (~ 12 kOe). The behaviour described by [4-6] is found in our simulations. When the 
field is decreased starting from the saturated up state (see fig. 2): down bubbles first appear, 
decreasing the field further, bubbles are elliptically deformed (stripe out), then coalesce to 
form a maze structure, the domain walls move progressively making up domain size decreases, 
and finally up bubbles appear. We may emphasize that the whole evolution is completely 
described using only the hamiltonian with no further hypothesis. 

Domains in a demagnetized state for different potential values are shown in fig. 3. The 
V+ part of the potential is arbitrarily fixed to unity since the calculation is performed under 
zero field. The demagnetized state is obtained starting from a full random initialization. This 
corresponds physically to demagnetize the film through annealing it at an infinite temperature, 
and to cool it down abruptly. We have tested that it is not necessary to use a complex 
procedure of simulated annealing to reach the ground state of the system. Our domain patterns 
in fig. 3 can be compared to the MFM images obtained by V. Gehanno et al. [3] , for different 
film thicknesses. The domain geometry evolves from a maze type structure like our fig. 3a to 
a fig. 3c type domain structure by decreasing the film thickness from 50 nm to 5 nm, i.e by 
decreasing the V- dipolar part of the potential. The pattern in fig. 3d fits well to SEMPA 
images obtained by R. Allenspach et al. [12] for ultra-thin cobalt films (~ 0.8 nm) on the 
(111) gold surface. Since do ~ 50 nm, the image size in fig. 3a-d is around 10 /nm x 10 /xm. 
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Fig. 4 - Schematic representation of the first, second and third layers of the i cell neighbours. The 
potential V is constant on each layer. 

Fig. 5 - C+ and C_ points are marked by respectively + and -. If a hyperplan (P) separating C+ 
from C- points is found, a vector W perpendicular to (P) verifies that the i^*-'^ projection on W is 
positive or negative depending belongs respectively to C+ or C- classes. 



which gives a domain size: 2-5 /im for fig. 3d to 50-100 nm for fig. 3a-c, in good agreement 
with experimental data [3,12]. 

We now solve the inverse problem: can we from a given calculated or experimental domain 
pattern, infer the unknown potential VI One needs to find the relation between the state of 
the cell i and envi. The i cell's environment can be described with the vector C whose 
the q components are c}q^ = X)pe<5<,(j) ^P' '^here the 5q{i) set states for the i cell's q layer of 
neighbours (see fig. 4). We also define a vector W with components Wq equal to the 
j cell belonging to 5q{i). We separate the C, vectors into two classes C+ and C_ depending 
on rji respectively equals -f 1 or -1. If one excepts the marginal case where = (which is 
only true for B ^ B^), the rule © reads: 

~C^"\w) + B > 

->,, ->( ~ (7) 

where (^C — J2q Cq^^g — J2k ^kVi+k is the euchdian scalar product. In other words, 

one has to find a hypcrplanc (P) separating C+ points from C_ points, thus a W vector 
perpendicular to (P) will be solution of eq. [71(the case S = is illustrated on fig. 5). 

We start with a randomly chosen vector W'^'^^ . For a given i cell, the C vector is 
calculated. If (w^°'^ , "C ^ -f i? and 77, have the same sign, another cell is tested. In the other 

case, is modified according to the rule 

ly (1) = W^") sign (^(W^"'> , c'(^) ) + s) ] T^^) (8) 

and so on until a vector iy(") is found, verifying 77^ = sign (J^W'^^^^ C '•^''^ + for every i 
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Fig. 6 - Potentials corresponding respectively to a,b,c,d parts of fig. 3. The guessed potentials are 
represented by points, the ideal ones by solid lines and the dashed lines are the smoothed guessed 
potentials. 



cell. Finaly, the components give the V potential we are looking for. 

In fact, this algorithm follows closely the Rosenblatt's perceptron NNT [13] in its simplest 
expression since only one neuron is needed. The N components Cg*^ are injected in the N cells 
input layer of the perceptron. The neuron answer a is sign (J^W''"^\ C ^'"'^ + i the Wq"^'' 
being the synaptic coupling strengths between the q cell of the input layer and the neuron. 
If the perceptron makes a mistake {rji ^ cr), the synaptic couplings are corrected following 
eq. IHl this equation defining the error-correction learning rule of the perceptron. Papert and 
Minsky [14] have shown that the perceptron NNT converges in a finite number of steps if 7 > 

and if the classification problem of the C vectors between the two classes C+ and C_ is a 
linearly separable one [15]. In our case, the separability condition is verified by construction 
and for 7 = 0.5, a correct potential is found after only one scan of a 250 x 250 cells domain 
pattern. 

The guessed potentials V jViX) from demagnetized domain structures of fig. 3 are shown 
in fig. 6. The shape of the potential (positive between first neighbours and negative for other 
neighbours) as well as i?_ and i?+ are estimated correctly. The determination of y_/V+ 
is less accurate because of statistical fluctuations. These latters have been reduced here by 
taking the mean value of the curve between neighbour points (dashed line). We have only 
used demagnetized states because it is often the only information that can be experimentally 
obtained, but one can use domain patterns corresponding to different values of B as an input 
and statistical fluctuations will be strongly reduced. Nevertheless, the above results show that 
even using a single image, and such a simple perceptron algorithm, it is possible to obtain 
probing results. 

In summary, magnetic bubble films can be described using a 2 D super-Ising hamiltonian. 
The potential allowing to fit an experimental or calculated domain pattern can be guessed by 
a perceptron neural network. The universality of the super-Ising model and of the bubble and 



Stephane Padovani: How to GUESS the inter magnetic bubble etc. 



7 



maze domains [16] suggests this approach could be exploited in numerous other problems as 
demixion processes [17], Turing mechanisms [18] or surface gases [19] to analyse experimental 
images. 

* * * 

I am very grateful to F. Bardou for fruitful discussions. 
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